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A  recent  multi-sensor  optical  turbulence  measurement  campaign  in  the  Haute-Provence  region 
of  France  provided  a  unique  comparison  of  optical  turbulence  sensors.  Thermosondes  from  both 
the  US  and  France  were  used,  measuring  temperature  structure  functions  for  30cm,  95cm  and  one 
meter.  In  addition,  turbulence  was  measured  using  a  Generalized  Scidar  (GS)  mounted  on  a  1.93m 
diameter  telescope  at  l'Observatoire  de  Haute-Provence  (OH.  The  GS  measures  optical  turbulence 
strength  and  altitude  from  the  scintillation  of  binary  stars  observed  over  periods  from  several 
minutes  to  several  hours.  The  GS  has  300m  vertical  resolution  from  25km  down  to  nearly  the  top 
of  the  telescope.  This  paper  describes  the  instruments  and  compares  their  performance.  While  the 
instruments  did  not  yield  identical  results,  agreement  was  reasonable. 


Nomenclature 


a 

= 

Constant 

r 

= 

Distance  for  structure  function  or  C**  (m) 

Bbs ** 

= 

Difference  in  autocorrelation  function 

r0 

= 

Fried’s  coherence  length 

b 

= 

Vertical  stability  of  atmosphere  (m'1) 

S 

= 

Autocorrelation  of  detector  impulse  response 

C 

= 

Autocorrelation  function 

Srr 

= 

Power  spectral  density  of  turbulence 

c„2 

= 

Index  of  refraction  structure  constant  (m’2  3) 

V 

= 

Velocity 

CT2 

= 

Temperature  structure  constant  (K.2tn'2/3) 

X 

= 

Distance  in  an  arbitrary  direction  (m) 

Dy(r) 

= 

Structure  function  of  variable  y 

a 

= 

Constant 

h 

= 

Altitude  above  ground  (m) 

X 

= 

Constant  for  balloon  motion  (s'1)  See  Eq.  (15) 

K 

= 

Turbulence  wave  number  (m_1) 

AH 

= 

Vertical  Resolution 

k 

= 

Wave  number  of  light  (mf1) 

Am 

= 

Difference  in  star  magnitude 

Z 

= 

Equivalent  width  of  autocorrelation  function 

e 

= 

Angular  separation  of  binary  star 

N 

= 

Noise 

en 

= 

Isoplanatic  angle 

P 

= 

Pressure  (Pa) 

A 

= 

Wavelength  (m) 

T 

= 

Absolute  temperature  (K) 

ct! 

= 

Rytov  varience 
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1.  Introduction 

It  is  a  well-known  hypothesis  that  turbulence  in  the  upper  atmosphere  is  often  the  result  of  gravity  (buoyancy)  waves 
propagating  into  the  upper  atmosphere,  growing  in  amplitude  with  altitude  until  they  become  unstable,  then  breaking  into 
turbulence.  Coulman  et  al.\  made  an  early  proposal  of  the  wave  -  turbulence  association  in  conjunction  with  the  reporting 
of  data  from  the  thermosonde,  a  balloon-borne  instrument  that  measures  turbulence,  which  is  described  below.  However, 
validation  of  the  wave  hypothesis  has  proven  to  be  difficult.  In  July  2002,  a  study  of  the  wave  -  turbulence  interaction  was 
performed  by  a  French  team,  lead  by  Coulman’s  colleague,  Jean  Vernin,  and  a  team  from  the  Air  Force  Research 
Laboratory,  Hanscom  AFB,  MA.  This  campaign  utilized  a  number  of  sensors  to  detect  gravity  wave  activity  and  optical 
turbulence.  Gravity  waves  were  sensed  by  the  radiosonde  instruments  on  the  thermosondes2.  Turbulence  was  sensed  by  the 
Generalized  Scidar  (GS)  and  the  thermosondes.  Reference  2  concentrates  on  the  extraction  of  gravity  wave  information 
from  the  balloon  data.  This  campaign  provided  a  unique  opportunity  to  compare  the  performance  of  the  optical  turbulence 
sensors,  which  is  the  focus  of  this  paper. 

A.  Thermosondes 

The  thermosonde3,4  has  been  used  in  optical  turbulence  research  since  1971.  It  is  attached  to  a  meteorological 
radiosonde  and  is  carried  into  the  atmosphere  by  a  balloon.  The  thermosondes  flown  by  Air  Force  Research  Laboratory 
personnel  use  two  very  fine  resistance  wire  probes  to  measure  the  temperature  difference  across  a  lm  horizontal  distance. 
Onboard  electronics  convert  the  temperature  difference  to  a  voltage  signal,  amplify  and  filter  the  signal,  then  perform  a  4  to 
8  second  running  average.  The  signal  is  then  transmitted  to  the  ground  station  along  with  the  standard  meteorological  data: 
temperature,  pressure,  humidity,  wind  speed  and  direction  sensed  by  the  radiosonde5,6. 

The  University  of  Nice  team  has  flown  thermosondes  with  a  common  origin  to  the  US  instruments,  but  which  evolved 
separately  through  an  early  Italian  design7.  The  French  thermosondes  use  4  wire  sensors  to  measure  two  separate  structure 
functions  with  lengths  of  0.3m  and  0.95m. 

1.  Thermosonde  Theory 

The  use  of  the  thermosonde  to  estimate  atmospheric  turbulence  is  based  on  the  statistical  theory  of  a  randomly  varying 
media.  If  two  temperature  sensors  pass  through  the  atmosphere  a  distance,  r,  apart  along  a  path,  x,  the  mean  value  of  the 
square  of  the  difference  is  defined  as  the  temperature  structure  function,  DT  (r ) 8: 

Dr(r)  =  ([r(x)-r(x+r)]2)  CD 

Structure  functions  of  different  lengths  generally  have  different  values.  The  theory  that  relates  the  different  sizes  comes 
from  Kolmogorov.  Kolmogorov  hypothesized  such  a  relationship  for  velocity  fluctuations  in  a  locally  homogeneous  and 
isotropic  random  field  for  scales  less  than  the  largest  eddies  known  as  the  outer  scale,  L0,  which  are  on  the  scale  of  the 
turbulent  kinetic  energy  source,  and  greater  than  the  smallest  eddies  known  as  the  inner  scale,  10,  which  are  the  scales  of 
molecular  energy  dissipation.  Kolmogorov  reasoned  that  energy  is  transferred  through  inertial  processes  from  the  outer 
scale  to  the  inner  scale,  the  range  between  them  is  known  as  the  inertial  range.  By  dimensional  arguments,  Kolmogorov 
reasoned  that,  at  equilibrium,  longitudinal  velocity  fluctuations  (along  the  direction  of  the  vector  r)  in  this  range  should 
satisfy  the  structure  function  relationship: 


Dv{r)  =  CyA  (2) 

where  the  constant  C2  is  known  as  the  velocity  structure  constant.  While  the  condition  of  inertial  range  “equilibrium 
turbulence”  defined  above  is  not  always  present  in  the  atmosphere,  functional  relationships  with  a  power  near  2/3  are 
surprisingly  prevalent  in  the  atmosphere.  The  2/3  power  in  the  structure  function  is  equivalent  to  the  one  dimensional 
power  spectrum  (for  positive  frequencies)  shown  in  Eq.  (3).  Since  the  power  of  the  spatial  wave  number  is  -5/3,  the 
Kolmogorov  profile  is  either  referred  to  as  the  -5/3  power  law  model  or  the  2/3  power  law  model. 


Srr{K)  =  0.25Cv2iT^  (3) 
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Velocity  turbulence  in  the  presence  of  temperature  gradients  causes  temperature  fluctuations.  Corrsin,  Yaglom,  and 
Obukhov9  independently  extended  Kolmogorov  theory  to  passive  scalars,  and  showed  that  the  structure  function  of 
conservative  passive  scalars  should  also  follow  the  2/3-power  law.  Therefore,  the  temperature  structure  function  is 
expressed  in  terms  of  temperature  structure  constant,  CT. 

Dr{r)  =  Clr%  (4) 

Cj.  is  not  sufficient  for  prediction  of  turbulence  disturbances  to  the  propagation  of  electromagnetic  radiation.  These 
disturbances  are  based  on  the  structure  constant  of  the  index  of  refraction  of  the  air,  Cn2.  Index  of  refraction  is  related  to  air 
density  through  the  Gladstone  Dale  Constant.  Turbulence  in  the  free  atmosphere  is  slow  enough  that  variations  in  density 
are  essentially  due  to  the  turbulent  fluctuations  of  temperature,  since  pressure  disturbances  dissipate  at  the  speed  of  sound. 
For  microwave  radiation,  the  variations  in  humidity  must  also  be  considered.  Conversion  of  the  temperature  structure 
constant  to  the  refractive  index  structure  constant,  Cn2,  depends  on  local  pressure  and  mean  temperature  and  the  wavelength 
of  the  radiation  that  is  being  propagated.  For  radiation  near  the  visible  spectrum,  and  when  the  moisture  contribution  can 
safely  be  ignored  (normal  except  right  over  a  body  of  water)  it  is  customary  to  use  the  formula10*: 

Cl  IT2)*  C2t  (5) 

where  ad(k)  =  79xlO‘8K/Pa  for  visible  and  near-infrared  wavelengths. 

The  theory  described  above  implies  a  long  run  through  a  stationary,  randomly  varying  temperature  field  in  an 
equilibrium  cascade.  For  the  thermosonde,  the  conditions  are  varying  as  we  ascend  at  anywhere  from  3  to  8  m/s  through 
what  may  or  may  not  be  stationary,  homogeneous,  equilibrium  conditions.  Therefore,  we  refer  to  our  results  as  an  estimate 
of  the  structure  constants. 

2.  Thermosonde  Hardware 

There  are  differences  in  the  way  the  U.S.  and  French  apply  Kolmogorov  theory  to  the  design  of  the  thermosonde.  The 
U.S.  thermosondes  have  sensors  1  meter  apart,  as  shown  in  Figure  1,  measuring  the  one-meter  structure  function.  The 
French  thermosondes  uses  four  probes  to  estimate  the  0.3m  and  0.95m  structure  functions,  as  shown  in  Figure  2.  They 
sometimes  use  other  lengths.  Eq.  4  provides  a  relationship  among  the  different  sized  structure  functions,  and  used  to 
compute  the  temperature  structure  constant,  Cj,  which  should  be  independent  of  distance.  The  U.S.  measurement  across  1 
meter  simplifies  the  process  since  the  structure  constant  is  equal  to  the  1-m  structure  function.  The  French  team  uses  the 
two  structure  functions  to  produce  two  estimates  of  the  structure  constant. 


Figure  1.  A  U.S.  thermosonde,  with  a  1-m  distance 
between  the  two  sensors  on  the  boom.  Radiosonde  is 
taped  in  front,  structure  function  electronics, 
batteries,  and  spare  channel  interface  are  behind. 


a  0.95m  distance  between  pairs  of  temperature 
sensors.  Radiosonde  electronics  will  be  infront  of 
container,  other  electronics  behind. 


*  In  this  reference,  one  of  the  square  symbols  (2)  was  inadvertently  omitted.  Please  correct  it! 
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A  second  design  difference  is  the  configuration  of  the  resistance  wires  used  to  sense  temperature.  The  resistance  wire 
for  the  U.S.  instrument  is  tungsten  with  a  diameter  is  3.45pm.  The  length  of  the  U.S.  wire  sensor  is  less  than  5mm, 
stretched  between  two  conducting  pins,  as  shown  in  Figure  3.  The  nominal  resistance  is  27f2,  and  the  temperature 
coefficient  of  resistivity  of  0.00375  K'1.  The  current  through  the  wires  is  approximately  0.1  A.  The  wires  are  approximately 
perpendicular  to  the  airflow  as  the  balloon  lifts  the  instrument  through  the  atmosphere;  similar  to  typical  “hot-wire”  and 
“cold-wire”  probes  in  wind  tunnels.  U.S.  sondes  are  typically  launched  with  a  target  ascent  speed  of  5  to  7m/s. 

In  the  French  design,  the  tungsten  wire  is  5pm  diameter  and  60mm  long.  It  is  threaded  between  7  looped  holders  to 
form  a  conically  shaped  crown  on  the  top  of  the  vertical  holder,  as  shown  in  Figure  4.  The  geometry  of  the  design  is 
intended  to  provide  a  sample  of  approximately  a  cubic  centimeter  of  air.  The  nominal  resistance  is  220F2  at  20°C  with  a 
resistivity  coefficient  of  0.0035  K  !.  The  current  through  each  wire  is  1mA.  The  target  ascent  speed  is  4m/s. 

In  both  designs  a  pair  of  sensor  wires  forms  two  legs  of  a  Wheatstone  bridge  that  generates  a  rapidly  varying  voltage 
based  on  the  difference  between  the  resistances  of  the  wires.  The  signals  are  amplified,  filtered,  and  then  averaged.  The 
U.S.  design  uses  a  circuit  that  performs  a  running  average  with  a  time-constant  of  about  4  to  8  seconds.  The  French  design 
uses  a  circuit  with  a  time  constant  for  averaging  of  about  a  1.25sl.  The  French  design  has  a  duplicate  circuit  on  the  board 
for  the  second  structure  function. 

Both  teams  currently  use  a  Vaisala  RS-80  Radiosonde  with  spare  channel  capability  to  send  the  running  average  of  the 
structure  function  and  some  ancillary  information  down  to  the  ground  station  with  the  typical  atmospheric  parameters 
sensed  by  weather  balloon  instruments:  ambient  temperature,  pressure,  relative  humidity,  and  wind.  Altitude  is  determined 
by  the  hydrostatic  equation,  corrected  for  humidity. 

There  is  an  interesting  difference  in  design  philosophy  between  the  two  thermosondes.  In  the  case  of  the  U.S.  design, 
the  field  operators  are  laboratory  employees,  often  engineers  and  technicians.  Consequently,  the  field  operations  have  been 
optimized  for  an  expert  staff  in  the  field  for  final  calibration  and  assembly  of  each  package.  In  the  case  of  the  French 
design,  the  operators  are  often  graduate  students,  and  the  packages  have  been  engineered  for  ease  of  field  operation. 
Consequently  much  of  the  calibration  is  performed  back  at  the  University  of  Nice  laboratory.  Probes  are  already  paired  to 
insure  nearly  identical  resistance,  color-coded  and  packaged  for  each  thermosonde  in  groups  of  4.  In  the  field,  the  probes 
are  fastened  to  the  thermosonde  with  simple  connectors  to  the  color-coded  receptors.  The  shielded  electronics  are  slipped 
into  a  cardboard  container  which  folds  into  an  aerodynamic  shape  ready  for  flight. 


Figure  3.  Sensor  wire  end  of  U.S.  probe.  Resistance 
wire  is  strung  between  the  two  electrodes. 


Figure  4.  Sensor  wire  end  of  French  probe. 
Resistance  wire  is  strung  from  a  top  loop  to  a  bottom 
loop  and  back  to  the  top  forming  a  conically  shaped 
crown. 
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B.  Generalized  Scidar 

A  scidar  is  an  optical  method  of  measuring  optical  turbulence".  Short  exposure  time  images  such  as  shown  in  Figure  5 
are  analyzed  to  determine  the  profile  of  optical  turbulence  in  the  beam  path.  In  the  classical  scidar,  images  from  the 
telescope  pupil  are  analyzed,  which  makes  it  insensitive  to  turbulence  close  to  the  ground  because  the  scintillation  variance 

is  proportional  to  ,  where  h  is  the  altitude  above  the  ground  of  the  turbulent  layer  (acting  as  a  phase  screen)12.  In  the 
GS  the  plane  of  the  detector  is  made  the  conjugate  to  a  plane  (analysis  plane)  at  a  distance  h  ,  on  the  order  of  a  few 

kilometers,  below  the  telescope  pupil  {hgs  >  0) .  With  a  GS,  the  turbulence  at  ground  level  (including  the  telescope  dome) 
becomes  detectable  because  the  distance  relevant  for  scintillation  produced  by  a  turbulent  layer  at  an  altitude  h  is  now 
|  h  —  hgs  |  above  the  analysis  plane. 

The  autocorrelation  of  the  scintillation  produced  by  a  turbulent  layer  consists  of  three  components:  one  is  centered  at 
the  origin,  and  the  two  others  are  separated  by  OH  and  —OH ,  respectively,  where  0  is  the  angular  separation  of  the 
double  star,  and  H  =  |/t  —  hgs  j ,  which  is  equal  to  h  —  hgs  for  the  case  of  interest,  since  hgs  <  0 .  As  the  different  turbulent 

layers  are  statistically  independent,  the  contribution  of  each  one  is  added,  and  the  total  theoretical  autocorrelation  function 
can  be  written  as: 


C“  (r)  =[dh  C2  ( h )  {«C(r,  h)  +  b  [C  (r  -  Oh,  h)  +  C  (r  +  Oh,  ft)]} 


(6) 


where  C  (r,  //)  represents  the  autocorrelation  of  the  scintillation  of  a  single  star  produced  by  a  layer  of  unit  C2 ,  at  an 
altitude  h .  The  factors  a  and  b  are  given  by 


a  =  - 


1  +  a 


b  = 


a 


(1  +  tz)2’  (1  +  a)2’ 


a  =  10 


,-0.4Am 


(7) 


where  Am  is  the  magnitude  difference  of  the  double  star. 

Eq.  (6)  shows  that  all  the  information  needed  to  retrieve  C 2  {hi)  is  contained  in  a  radial  section  of  C  (r)  along  the 

double  star  separation.  Furthermore,  we  wish  to  eliminate  the  central  peak  where  the  contribution  of  each  layer  is 
undistinguishable  from  that  of  the  others,  as  they  are  added,  and  it  contains  the  uncorrelated  noise.  For  this  reason,  the 
difference  of  the  sections  of  the  measured  autocorrelation  function  parallel  and  perpendicular  to  the  stars  separation, 
C|p  and  C;  respectively,  are  calculated.  This  difference  can  be  written13: 


luw-cr-cr 

=  J  dhK(x,h)C2„(h-hg5)  +  N(x) 


where  N(x) is  the  noise,  and  the  kernel  K (x, h)  is  a  radial  section  of  C(y  —  0H,H^*S(j)  along  the  direction 
r  ||  <9 ,  where  S'(r)  is  the  autocorrelation  of  the  impulse  response  of  the  detector.  By  measuring  Bgs„  {x) ,  calculating 
theoretically  K  (x,  ti)  and  knowing  hgs  as  well  as  an  estimate  of  the  noise  variance,  Eq.  (8)  (of  Fredholm  type)  is  inverted 
to  retrieve  Cfj  (h)  using  a  maximum  entropy  algorithm.  The  inversion  is  simplified  by  the  diagonal  elements  in 
K  (x,  h)  along  the  line  x  =  0  [h  -  hgs ) . 
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Important  features  of  the  GS  should  be  noted:  The  space  between  the  telescope  pupil  and  the  analysis  plane  is 
turbulence-free  because  it  is  virtual.  The  C2  ( h )  measurements  are  independent  of  static  optical  aberrations,  as  they  are 

cancelled  out  in  the  data  analysis.  The  vertical  resolution  of  the  retrieved  C„  (//)  profiles  is  given  by14 


A H  = 


z(M<.,|) 

e 


(9) 


where  X|J/t  —  /2gJ,|j  is  the  equivalent  width  of  the  scintillation  spatial  autocorrelation  function,  related  to  the  first  Fresnel 
zone,  and  A  is  the  wavelength.  As  an  example,  the  resolution  achieved  at  ground  level  (/t  =  0)when  observing  the 

double  star  95  Herculis  (  0  =6.2  arcsec),  with  the  analysis  plane  at  [hgs  =  -4km} ,  is  A  FI  =  740m.  But  the  maximum 

entropy  method  increases  the  vertical  resolution  by  a  factor  of  about  2,  which  leads  to  a  final  vertical  resolution  of  about 
300m. 

A  typical  C2n  plot,  shown  in  Fig.  6,  is  obtained  using  binary  star  y  Del  with  angular  separation  of  9  -  9.3 

arcsec.  Many  layers  are  detected  at  1km,  just  above  ground  level,  at  4km,  and  then  4  layers  are  visible  between  10  and 
15km.  The  vertical  resolution  is,  as  expected,  about  300m.  The  temporal  resolution  is  about  21s. 

The  GS  performance  depends  on  the  size  of  the  telescope  mirror.  The  separation  distance  of  the  two  star  images  is 
equal  to  the  angular  separation  of  the  stars  times  the  altitude  of  the  turbulence  layer.  The  factors  that  limit  the  maximum 
angular  separation  are:  the  images  migrate  on  the  image  based  on  the  wind  velocity  of  the  turbulent  layer,  so  the  distance 
between  the  outer  images  must  be  limited  to  a  fraction,  /,  of  the  mirror  diameter;  and  it  is  desired  to  measure  turbulence 
layers  up  to  25km.  The  maximum  separation  is  6mm  —  f  D I  hmm  where  D  is  the  mirror  diameter.  The  new  very  large 

telescopes  permit  the  use  of  much  greater  star  separations,  opening  up  the  opportunity  for  star  pairs  that  are  not  necessarily 
binary  stars,  permitting  longer  coverage  times  at  higher  zenith  angles. 


1.6-C110  19-  1-9S  5:10  Til  EtaO.  daubb  At  m  20  m 


-JO  0  SO 

V„  m/ti 


Figure  5.  Cross-correlation  of  scintillation  images 
separated  by  20ms,  obtained  with  the  double  star 
KPupis,  of  angular  separation  of  9.9  acrsec.  The 
coordinates  of  the  central  peak  of  each  of  the  3 
triplets  give  the  wind  velocity  of  the  corresponding 
layer. 


G-SCOW  193-GHP-  1 8/  7/  2  0:47->  2:  9  U T  <SEflND>=0.85  mtlM  N 


Figure  6.  Graphical  display  of  GS  results  showing  a 
time  series  of  optical  turbulence  profiles  up  to  25km 
above  OHP  as  a  function  of  time  for  a  period  from 
before  1:00  to  after  2:00UTC  on  18  July  2002.  The 
color  scale  depicts  C„2  levels  from  ~lx!0'18  (white)  to 
lxlO16  m'J/3  (dark  blue). 
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II.  The  Experiments 

The  campaign  started  the  evening  of  17  July  2002,  and  lasted  until  the  morning  of  24  July  2002.  Data  were  taken  from 
two  separate  locations.  GS  optical  data  and  thermosonde  data  were  acquired  at  I’Observatoire  Haute  Provence  (OHP) 
(latitude  43 .9N  and  longitude  5.7 IE  and  altitude  630m  MSL).  Thermosonde  data  were  also  acquired  at  the  private 
SIRENE  Observatory  located  at  Lagarde  d’Apt  (designated  in  this  paper  as  “APT”)  (latitude  44. ON  and  longitude  5.5E  and 
altitude  1,100m  MSL). 

The  Generalized  Scidar  (GS)  was  mounted  on  the  1.93m  diameter  telescope  at  OHP.  The  GS  provides  continuous 
nighttime  profiles  of  turbulence  at  300m-resolution  from  telescope  altitude  up  to  25km.  Thermosondes  were  launched  at 
the  OHP  site  and  at  the  APT  site,  which  is  20km  W  by  NW  of  OHP.  Gravity  wave  activity  was  deduced  from  temperature, 
horizontal  wind  velocity  and  balloon  ascent  rate  measured  by  radiosondes,  which  are  part  of  the  thermosonde  system. 
Over  the  course  of  the  week,  the  US  team  launched  8  sondes  from  APT,  the  French  team  launched  23  thermosondes  from 
APT  and  19  from  OHP,  and  acquired  4317  GS  files  of  approximately  20s  duration  each. 

On  the  night  of  19-20  July,  two  flights  were  made  with  a  U.S.  and  French  thermosonde  attached  to  the  same  balloon. 
In  one,  the  thermosondes  were  secured  next  to  each  other  with  the  booms  at  equal  altitude  50m  below  the  balloon.  In  the 
second  flight  the  French  thermosonde  was  located  50m  below  the  balloon  and  the  U.S.  thermosonde  was  located  50m 
below  the  French  instrument. 


III.  Analysis  and  Results 


A.  Thermosondes 

I.  US  and  French  Thermosonde  Comparison 

The  first  dual  thermosonde  flight,  with  the  US  and  French  instruments  right  next  to  each  other  50m  below  the  balloon, 
provided  some  very  worthwhile  data  for  comparison;.  The  second  flight  with  the  US  instrument  50m  below  the  French 
instrument  did  not  work  as  well.  There  appeared  to  be  interference  affecting  the  French  instrument,  especially  at  the  higher 
altitudes.  Therefore,  only  the  first  flight  will  be  used  for  comparison.  Because  of  the  differences  in  the  signal  averaging 
time  between  the  US  and  the  French  instruments,  the  actual  time  series  of  profile  of  data  appear  quite  different.  Therefore, 
for  purposes  of  comparison,  each  Cn2  profile  has  been  averaged  over  300m  bins.  The  results  for  the  Cn2  from  the  two 
French  structure  functions  (30cm  and  95cm)  and  the  single  lm  structure  constant  from  the  US  instrument  are  shown  in 


Figure  7.  With  a  few  regions  of  disagreements,  the 
comparison  is  quite  encouraging.  The  agreement 
between  the  US  (lm)  and  the  FR  95cm  is  very  close, 
which  might  be  expected  since  these  probes  were  very 
close  together.  The  C„2  from  the  30cm  data  is  also 
close  to  the  other  data,  except  for  the  region  between  8 
and  12km,  where  these  data  exceed  the  results  from  the 
larger  structure  functions  by  quite  a  bit.  This  could  be 
the  result  of  some  very  non-Kolmogorov  behavior,  say 
a  region  of  very  low  turbulence  that  still  has  a  large 
amount  of  very  small-scale  turbulence.  However,  it  is 
not  a  good  idea  to  put  too  much  credence  on  the  results 
from  a  single  flight.  It  is  possible  that  the  electronics 
for  that  structure  function  experienced  a  high  noise 
condition  for  that  altitude  range.  There  are  a  few  places 
where  the  US  instrument  is  showing  some  regions  of 
higher  turbulence  that  the  FR  instrument  does  not 
agree:  for  instance  near  18km  and  near  20km. 

The  impact  of  turbulence  on  light  passing  through 
the  atmosphere  is  a  function  of  integrals  through  the 
region.  Some  optical  performance  integrals  are:  The 

Rytov  variance,  cr  a  measure  of  the  amount  of 

scintillation  of  the  light;  Fried’s  Coherence  length,  r0 
an  indication  of  the  largest  uncorrected  telescope  lens 
that  can  be  used  effectively;  the  isoplanatic  angle,  0{} ,  a 


Figure  7.  Comparison  of  C„2  computed  from  French 
95cm,  30cm,  and  US  lm  structure  functions.  Data  have 
been  averaged  over  300m  bins.  Balloons  launched  from 
Apt  on  19  July  2002. 
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measure  of  the  maximum  angle  of  that  one  could  expect  coherence  between  an  incoming  beam  from  a  source  and  an 
outgoing  beam.  The  latter  would  apply  to  the  ability  to  use  adaptive  optics  to  correct  observed  phase  distortions.  For  light 
coming  down  through  the  atmosphere  from  a  source  at  great  distance,  say  a  star,  the  equations  for  these  three  performance 
parameters  are: 


0-1=0.56^  dTj 


ro  = 


0.423  k 


iCl{v)dri 


-X 


#0  = 


2.9U2fc„2(/7)7K  dr, 


-X 


(10) 

(11) 

(12) 


where  k  is  the  wave  number  of  the  light  (  k  =  2nj X ),  Tj  is  the  length  along  the  beam  path,  and  L  is  the  total  beam  path. 
Each  integral  is  terminated  at  an  altitude  where  the  C„2  is  so  small  that  it  no  longer  contributes  to  the  integral. 

The  three  performance  parameters  were  evaluated  for  the  three  Cn2  profiles.  There  were  some  rather  large  differences 
in  the  three  C„  results  near  the  ground,  where  turbulence  is  very  high,  but  high  reading  can  come  from  instrument,  so  the 
performance  integrals  were  evaluated  from  2km  up  to 

Table  1.  Comparison  of  optical  performance  parameters 
from  the  three  sources  of  Cn2,  the  US  lm  structure 
function,  and  the  French  30cm  and  95cm  structure 
functions. 


21km  -  the  maximum  altitude  of  data  from  the  US 
instrument.  The  results  are  shown  in  Table  1.  The 
agreement  between  the  US  and  the  French  instruments 
are  reasonably  close,  as  expected  from  the  C„2  profiles. 
The  FR  30cm  result  is  worse  for  some  quantities,  but 
not  bad  in  others.  Since  each  performance  parameter 
has  different  weighting  with  distance,  the  results 
respond  accordingly.  Whereas  r0  is  equally  weighted 


o 

with  altitude,  the  weighting  for 
altitude,  and  the  weighting  for 


O0  decreases 
_ 2  . 

cr  1  increases 


with 

with 


Cnl  Source 

r0  (cm) 

0o{jUmd) 

(n.d.) 

US  lm 

25.31 

11.58 

0.019 

FR  30cm 

20.80 

10.95 

0.023 

FR  95cm 

25.00 

13.59 

0.016 

altitude  until  the  integration  stops  due  to  insignificant 
C  2 

'-'ll  • 

2.  Kolmogorov  Hypothesis 

As  mentioned  above,  the  French  use  the 
Kolmogorov  hypothesis  to  compute  two  estimates  of 
CT2  from  the  two  structure  functions, 

Dt  (.3)  and  DT  (-95) ,  using  Eq.  (4).  Many  groups 

have  studied  the  applicability  of  Kolomogorov’s 
hypothesis  for  turbulence  in  the  atmosphere.  Two 
groups  have  compared  many  structure  functions 
obtained  from  ascending15  or  descending16,17  balloons. 
The  French  thermosondes  provides  another  opportunity 
to  test  the  applicability  of  the  hypotheses  to 
atmospheric  temperature  turbulence.  A  cross-plot  of 

Cn2  from  Dr  (.3) ,  designated  as  C„2-30,  to  Cn2  from 
Dt  (-95)  Cn2-95  for  each  time  is  shown  in  Figure  8. 

The  straight  line  is  the  identity  line.  A  first  impression 
is  that  the  points  generally  follow  the  Kolmogorov 
trend,  with  a  spread  in  the  data  of  about  an  order  of 
magnitude  about  identity.  It  also  appears  that  the 


Figure  8.  Cross-plot  of  C„2  from  the  30cm  structure 
function  to  C„2  from  95cm  structure  function  from  each 
data  timo  linp. 
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lowest  values  of  Cn2  might  have  higher  values  from  the 
30cm  function.  These  C„~  values  could  be  analyzed 
further;  however,  this  could  be  misleading.  Recalling 
the  strong  atmospheric  weighting  when  converting 
from  CT2  to  Cn2  from  Eq.  (5),  it  is  clear  that  the  altitude 
plays  a  very  strong  roll  in  the  lowering  Cn2  as  altitude 
increases.  Therefore,  there  is  strong  altitude 
dependence  in  the  placement  on  the  abscissa  of  the  plot. 

An  alternate  approach  is  to  compare  temperature 
structure  constants,  CT2  from  each  probe  pair.  The 
result  is  shown  in  Figure  8.  Again,  there  is  an 
impression  of  general  agreement  with  the  identity  line, 
and,  perhaps,  an  even  stronger  impression  that  at  low 
values  of  CT2,  the  values  for  CT2-30  are  larger  than  the 
corresponding  values  from  the  CT2-95.  Next,  the  ratio 
of  the  two  values  of  CT2  for  each  time  step  (CT2-30  / 
Ct2-95)  is  compared  to  CT2-95  in  Figure  9.  Here,  it  is 
quite  obvious  that  CT2-30  exceed  CT2-95  at  lower 
magnitudes,  but  that  the  ratio  decreases  as  the  CT2 
values  increase,  crossing  unity  near  the  middle  of  the 
population,  but  continue  to  decrease  as  Ct2  increases. 

Exploring  this  further,  the  data  were  separated  into 
bins  of  half  an  order  of  magnitude.  Next  the  average 
value  of  the  ratio  was  computed  for  each  bin,  and  the 
percent  of  the  total  population  of  3870  pairs  of  Cn2 
values  was  computed  for  each  bin.  These  results  are 
shown  in  Table  2  and  Figure  10.  This  appears  to  be  an 
almost  log  normal  population  centered  in  the  bin: 
lxl 0  S<  Cn2  <  lxlO'45  =  3.6xl0's.  It  is  clear  that  the 
average  ratio  is  quite  high  at  the  smallest  values  of  Cn2 
but  descends  to  near  unity  by  10-4  to  10-3.5,  remains 
near  there  as  C„2  increases,  until  the  final  bin,  where  it 
reduces  to  0.58.  There  are  only  5  points  in  this  bin,  so 
it  is  hardly  a  statistical  sample,  and  is  probably  not 
significant.  One  interpretation  of  these  results  is  that  for 
the  lowest  levels  of  turbulence,  the  larger  scales  have 
significantly  diminished  while  there  remains 
measurable  activity  in  the  very  small  scales.  At  this 
point,  this  is  mere  speculation. 

Table  2.  The  ratio  of  CT2  based  on  the  30cm  structure 
function  to  that  of  the  95cm  structure  function  for  bins 
of  half  an  order  of  magnitude  each  from  10'6  5  to  10"2. 
Also  shown  is  the  percent  of  total  population  in  each 
bin. 


Bin  Min 

(Iogio) 

Bin  Max 
(logio) 

Avg.  Ratio 
.3m  to  .95m 

Population 

(%) 

-6.5 

-6 

4.71 

2.2 

-6 

-5.5 

3.22 

10.9 

-5.5 

-5 

2.31 

21.3 

-5 

-4.5 

1.8 

27.6 

-4.5 

-4 

1.46 

20.6 

-4 

-3.5 

1.18 

11.7 

-3.5 

-3 

0.903 

4.4 

-3 

-2.5 

0.899 

1.08 

-2.5 

-2 

0.58 

0.13 

Figure  8.  Comparison  of  CT2  from  the  30cm  structure 
function  to  Ct2  from  95cm  structure  function  from  each 
data  time  line. 


1.E-07  1.E-06  1.E-05  1  .E-04  1  E-03  1.E-02 

Cr:  from  95cm  St.  Fen. 


Figure  9.  Comparison  of  the  ratio  of  CT2  from  the  30cm 
structure  function  to  CT2  from  95cm  structure  function 
from  each  data  time  line  vs.  CT2  from  the  larger  function. 


Figure  10.  Average  of  ratio  of  CT2  from  the  smaller 
structure  function  to  the  larger  vs.  logi0  of  CT2  binned 
by  half  magnitudes  (diamonds).  Unity  is  shown  as  the 
solid  line  near  the  bottom.  Also  shown  is  the 
population  distribution  in  the  bins  as  percent. 
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By  analyzing  Ct2  instead  of  Cn2,  does  not  imply  that 
the  range  of  CT2  is  independent  of  altitude.  A  plot  of  all 
the  CY  values. from  the  flight  we  have  been  analyzing, 
designated  VOL267,  is  shown  in  Figure  11.  The  values 
from  the  larger  structure  function  are  shown  in  dark 
diamonds,  and  those  from  the  smaller  are  shown  in  light 
circles.  The  minimum  CT2  values  recorded  are  typically 
the  noise  floor  of  the  instrument  that  does  not  directly 
depend  on  altitude,  but  the  electronic  components  of  the 
instrument,  which  may  be  dependent  on  altitude  related 
quantities  such  as  the  temperature  of  the  electronics 
compartment.  In  Figure  11,  it  appears  that  the  noise  floor 
of  the  larger  structure  function  is  lower  than  for  the 
smaller  structure  function.  The  maximum  values  of  Ct2 
are  generally  a  function  of  the  slope  of  the  potential 
temperature  profile  and  the  strength  of  the  turbulence. 

The  higher  values  at  low  altitudes  are  indicative  of  the 
high  turbulence  levels  near  the  surface,  and  the  generally 
higher  values  in  the  stratosphere  relate  to  the  increase  in  Figure  11.  Profiles  of  CT2  from  the  smaller  structure 
the  gradient  of  the  potential  temperature.  function  (small  lighter  circles)  and  from  the  larger 

These  results  are  from  one  balloon  ascent  only,  and  (darker  diamonds)  for  flight  VOL267. 
are  not  meant  to  imply  a  complete  analysis.  Rather  this 

has  been  presented  to  give  the  reader  an  appreciation  for  the  wealth  of  information  available  from  the  collection  of  more 
than  500  flights  of  the  French  thermosondes. 

B.  Scidar  and  Thermosondes 

A  comparison  of  a  time  series  of  optical  turbulence  for  the  GS  is  shown  next  to  a  Cn2  profile  from  a  thermosonde 
launched  early  in  the  period  shown  in  the  GS  plot.  The  output  is  in  fair  visual  agreement,  particularly  in  the  lower 


1.E-07  I.E-06 


t.M5  I.E-Qt  1.E-03 

CT2  SKV'4’') 


I.E-02 


•  e.-o  «a  >e-a  >  r  «  Cra>  a  >e.-o  >a 

Figure  12.  A  comparison  of  C„2  as  measured  by  the  GS  (left)  and  from  a  thermosonde  (smoothed  with 
a  300m  window).  The  GS  data  is  a  time  series  over  a  period  of  an  hour  and  20  minutes  on  the  evening 
of  18  July  2002.  The  thermosonde  was  launched  during  the  period  of  the  GS  data. 
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Trajectories  for  Evening  of  IS  July 
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Figure  13.  The  trajectories  of  5  thermosondes 
launched  on  the  evening  of  18  July. 


Figure  13.  A  comparison  of  the  C„2  measured  by 
the  GS  (solid  line)  located  at  the  altitude  of  the 
lower  dashed  line  and  2  thermosondes  launched 
from  the  upper  dashed  line14.  One  thermosonde  is 
shown  by  circles,  the  other  by  squares. 


atmosphere.  Both  plots  show  high  values  near  the  surface,  a  higher  band  from  2  to  3km,  a  band  between  4  and  5km,  and 
the  heavy  band  near  the  tropopause  between  troposphere  and  stratosphere  know  for  high  levels  of  turbulence. 

Above  the  tropopause,  the  comparison  breaks  down.  This  is  at  least  in  part  due  to  the  fact  that  the  thermosondes  are 
drifting  away  from  their  launch  sites  with  the  horizontal  velocity  of  the  prevailing  wind.  The  trajectories  of  5  balloons 
launched  that  night  are  shown  in  Figure  13.  The  telescope  is  located  at  the  origin  of  the  plot.  Since  the  GS  is  usually 
observing  double  stars  at  a  zenith  angle  of  less  than  45°,  the  outer  bounds  of  the  GS  data  would  be  a  circle  of  25km  radius 
drawn  around  the  origin. 

A  quantitative  comparison  of  the  measurements  of  the  GS  to  those  of  the  thermosonde  is  shown  in  Figure  14  from 
reference  14.  The  figure  compares  the  GS  located  at  the  lower  of  the  two  dashed  lines  to  two  balloons  from  a  campaign  in 
Chile  in  1998.  The  comparison  is  generally  good,  with  some  discrepancies  at  the  higher  altitudes,  probably  due  to  the  drift 
of  the  balloon  described  above. 


IV.  Conclusion 

The  2002  Campaign  in  the  Flaute-Provence  region  of  France  provided  an  excellent  opportunity  to  compare  turbulence 
results  from  three  methods  of  measuring  atmospheric  optical  turbulence:  the  generalized  scidar  (GS),  and  two  versions  of 
the  thermosonde,  the  design  used  by  the  French  team  and  that  used  by  the  US  team.  In  one  flight,  the  two  different 
thermosondes  were  lashed  together  and  flew  50m  beneath  the  ascending  balloon.  The  results  of  Cn2  as  determined  by  the 
lm  structure  function  of  the  US  instrument  closely  matched  the  Cn2  data  from  the  0.95m  structure  function  of  the  French 
instrument.  The  data  from  the  0.3m  structure  function  of  the  French  instrument  generally  agreed  quite  well  also,  except  for 
a  portion  of  the  atmosphere  where  it  showed  considerably  higher  turbulence  than  the  other  two  readings  for  an  unknown 
reason.  Optical  performance  parameters  from  the  three  structure  functions  generally  agreed  quite  well,  especially  Fried’s 
coherence  length  for  the  0.95m  and  lm  structure  functions  which  only  differed  by  1.25%. 

The  presence  of  two  structure  functions  on  the  French  instrument  provides  an  opportunity  to  test  Kolomogorov’s 
hypothesis  of  the  relationship  among  structure  functions  of  different  sizes.  Only  data  from  one  thermosonde  was  analyzed. 
The  analysis  showed  fair  agreement  with  the  hypothesis,  but  showed  increasing  divergence  at  the  lower  values  of  structure 
function,  with  the  smaller  structure  function  yielding  higher  than  the  values  predicted  by  the  hypothesis.  This  result  from  a 
single  ascent  demonstrates  the  wealth  of  data  that  is  routinely  available  from  the  French  design  to  more  fully  explore  the 
applicability  of  the  Kolmogorov  hypothesis  to  various  atmospheric  conditions. 

The  agreement  of  the  GS  with  the  balloons  is  especially  good  at  the  lower  altitudes.  At  higher  altitudes  the  agreement 
is  not  as  good,  but  the  thermosondes  typically  drift  out  of  the  region  of  GS  viewing  and  are  probably  experiencing  different 
atmospheric  conditions.  The  GS  is  a  remarkable  instrument  which  provides  nearly  continuous  profiling  of  the  turbulence 
in  the  atmosphere,  except  when  blocked  by  cloud  cover  or  when  binary  stars  are  not  in  the  viewing  region.  Measurement 
by  GS  is  clearly  the  instrument  of  choice,  at  least  whenever  a  telescope  of  sufficient  diameter  is  available. 
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